EPJ manuscript No. 

(will be inserted by the editor) 



Eikonal representation in the momentum-transfer space 



o 
o 

> 
O 

o 

cn 

m 
> 

(N 

m 
o 

Oh: 

I 

Oh, 



X 



P.A.S. Carvallio\ A.F. Martini^, M.J. Menon^ 

^ Centro Universitario de Patos de Minas, 38702-054 Patos de Minas, MG, Brazil 

^ Institute de Fi'sica Gleb Wataghin Universidade Estadual de Campinas, 13083-970 Campinas, SP Brazil 



Received: date / Revised version: date 

Abstract. By means of empirical fits to the differential cross section data on pp and pp elastic scatter- 
ing, above 10 GeV (center-of-mass energy), we determine the eikonal in the momentum - transfer space 
[q^- space). We make use of a numerical method and a novel semi-analytical method, through which the 
uncertainties from the fit parameters can be propagated up to the eikonal in the q^- space. A systematic 
study of the effect of the experimental information at large values of the momentum transfer is developed 
and discussed in detail. We present statistical evidence that the imaginary part of the eikonal changes sign 
in the g^- space and that the position of the zero decreases as the energy increases; after the position of 
the zero, the eikonal presents a minimum and then goes to zero through negative values. We discuss the 
applicability of our results in the phenomenological context, outlining some connections with nonpertur- 
bative QCD. A short review and a critical discussion on the main results concerning "model-independent" 
analyses are also presented. 

PACS. 13.85.Dz Elastic scattering - 13.85.-t Hadron-induced high-energy interactions 



1 Introduction 

High-energy particle scattering is the main experimental 
tool in the investigation of the inner structure of mat- 
ter. For particle-particle and antiparticle-particle scatter- 
ing, the highest energies reached in accelerators concern 
proton-proton (pp) and antiproton-proton (pp) collisions 
and these hadronic processes are expected to be described 
by the Quantum Chromodynamics (QCD), the gauge in- 
variant quantum field theory of the strong interactions. 
Recently, Dokshitzer stated that "QCD nowadays has a 
split personality. It embodies hard and soft physics, both 
being hard subjects and the softer the harder" In fact, 
despite of all the success of QCD in the treatment of hard 
and semi-hard scattering processes (large and medium mo- 
mentum transfer, respectively), the increase of the cou- 
pling constant when going to the soft sector (small mo- 
mentum transfer) does not allow the use of the most im- 
portant tool in theoretical physics, namely the perturba- 
tive calculation. On the other hand, the nonperturbative 
QCD approach, characterized by the non-trivial investi- 
gation of the vacuum structure and intricate simulation 
techniques in the study of bound states, did not yet pro- 
vide any result for the scattering states, based exclusively 
on its own foundations. Soft scattering embodies diffrac- 
tion dissociation (single and double) and elastic scattering 
12] and the point is that, presently, we do not know how to 
calculate even the elastic scattering amplitude (the sim- 
plest soft process) in a purely nonperturbative context. 
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Certainly, there are QCD-based or QCD-inspircd models, 
but, anyway, they have only bases or inspiration in QCD 

m 

However, from the experimental point of view, the ex- 
pectations in the field are very high due to the new gen- 
eration of colliders. From experiments that are being con- 
ducted at the BNL Relativistic Heavy Ion Collider (RHIC) 
it is expected data on pp soft scattering at center-of mass 
energies ^/s: 50 - 500 GeV, and in the near future the 
CERN Large Hadron Collider (LHC) will provide data on 
soft pp scattering at 14 TeV. 

At this stage, beyond the efforts directed to pure the- 
oretical developments (nonperturbative QCD) and phe- 
nomenology, "empirical" analyses play an important role 
in the search for model-independent results that can con- 
tribute with the establishment of novel and useful the- 
oretical calculational schemes. In order to achieve that, 
any approach must be based on the General Principles 
and Theorems of the underlying field theory, namely Uni- 
tarity, Analyticity, Crossing and their well founded con- 
sequences. Among the wide variety of formalisms the 
eikonal representation is distinguished by its intrinsic con- 
nection with Unitarity 5 and its efhciency as a useful 
framework for phenomenological analyses. Eikonal models 
can be distinguished by the different forms of the eikonal 
in the momentum-transfer space (q^ - space); for our pur- 
poses, we list some important and well known results in 
Refs. l6p7t.8^9 ,.10..imT^ . In this context, analyses based 
on empirical fits to the physical quantities involved and 
aimed to extract the characteristics of the eikonal as func- 
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tion of the momentum transfer, energy and reaction play 
an essential role in the attainment of consistent results in 
a model-independent way. 

In a previous work we have investigated pp scatter- 
ing at the CERN Intersecting Storage Ring (ISR) energy 
region, namely y/s = 23.5, 30.7, 44.7, 52.8 and 62.5 GeV. 
By means of fits to the differential cross section data and 
through a novel "semi-analytical" method, which allowed 
the error propagation from the fitted parameters, we have 
determined the imaginary part of the eikonal in the mo- 
mentum transfer space and have found statistical evidence 
for the existence of eikonal zeros (change of sign) in the 
region of momentum transfer = 7 ± 2 GeV^. These 
results have been used and discussed in different contexts 



In the present work, we have extended our previous 
analysis in several ways: (1) we include the pp data at 
lower energies, namely y/s =13.8 GeV and 19.5 GeV and 
also pp data in the region 14 GeV < < 1.8 TeV; (2) 
we make use of a numerical calculation in order to check 
the results obtained with the "semi-analytical" method; 
(3) we present in detail the statistical regions of confi- 
dence and all the numerical information; (4) we discuss 
the applicability of the results in the phenomenological 
context and also some possible connections with nonper- 
turbative QCD results. In addition, we present a critical 
discussion on the different analyses aimed to extract the 
eikoiral from fits to the differential cross section data. Our 
main novel conclusions are that the position of the eikonal 
zero decreases as the energy increases and, after the zero, 
the eikonal has a minimum and then goes to zero through 
negative values. 

The manuscript is organized as follows. In Sect. 12 we 
treat the Eikonal representation, recalling the formulas 
connecting the physical quantities that characterizes the 
elastic scattering with the eikonal in the q^- space and 
stressing the importance of the "empirical" information. 
In Sect. El after introducting the experimental data to 
be analyzed, we shortly review some typical results and 
open problems associated with model-independent anal- 
yses and, based on these considerations, we present our 
strategies, mainly related to the definition of two ensem- 
bles of experimental information. In Sect.Qlwe display the 
results of the fits to the differential cross section data, dis- 
cussing in detail our choice for the parametrization, the fit 
procedure and the effect of the uncertainties. In Sect.[Slwe 
treat the determination of the eikonal in the momentum 
transfer space by means of both the numerical and semi- 
analytical methods. In Sect. we discuss the applicability 
of our results in a phenomenological context and possible 
connections with nonpertubative QCD results. The con- 
clusions and some final remarks are the contents of Sect. 

m 

2 Eikonal Representation 

In this Section we first recall the essential formulas con- 
necting the eikonal in the momentum-transfer space and 
the physical quantities to be investigated. Next, we discuss 



the importance of the model- independent analyses, in par- 
ticular some questions related with eikonal zero (change 
of sign in the - space). 

2.1 General formalism 

In the Eikonal Representation, the elastic scattering am- 
plitude, as function of the center-of-mass energy, y/s, and 
the four-momentum transfer squared, = —t, is ex- 
pressed by 

F{s,q)=t 6dWo(g6){l-e'*(^^^)}, (1) 
Jo 

where b is the impact parameter, Jq is the zero order Bessel 
function (azimuthal symmetry assumed) and xi^i b) is the 
eikonal function in the impact parameter space. Although 
this formula may be deduced, for example, from the par- 
tial wave solution of the Schroedinger equation in the high 
energy limit "W , it should be stressed its character of well 
founded mathematical representation, independent of any 
approximation. In fact, as demonstrated by M. Islam, the 
complex scattering amplitude F{s,q) has a well defined 
representation in the impact parameter space, valid for all 
energies and scattering angles ^H]- This representation is 
usually named Profile function and denoted r{s,b). Since 
the exponential function is an entire function of its argu- 
ment, we may always express the Profile function in terms 
of an eikonal function: 

r{s,b) = i-e:^^''''\ (2) 

Unitarity allows to connect elastic and inelastic pro- 
cesses, which, in the impact parameter space, may be ex- 
pressed by the formula 

2Rer(s,6) = |r(s,6)|2+G,„ei(s,6), (3) 

where Gmeiis, b) is the Inelastic Overlap Function, namely 
the probability for an inelastic event to take place at s and 
b. This probabilistic interpretation demands that 

Ginel{s,b) < 1, 

and, from Eqs. (2) and (3), it is automatically verified for 

Im x(s, b) > 0. (4) 

This is the main result that characterizes the Eikonal rep- 
resentation, that is, the automatic agreement with Uni- 
tarity, a principle that, certainly, can never be violated. 

The Eikonal approach has a long history of success- 
ful results that allowed important developments. For ex- 
ample, from the earlies Chou-Yang iSj and Glauber PE] 
models, passing through the interesting geometrical for- 
malism by Bourrely, Soffer and Wu jl()| and, more recently. 
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the suggestive connections with quarks and gluons in the 
QCD-based models |11II12| . the Eikonal picture seems to 
represent a nearly natural framework and suitable phc- 
nomenological laboratory. 

All these eikonal models (and many others) may be dis- 
tinguished or classified according to different forms for the 
eikonal function in the momentum transfer space, which 
we shall represent by xi^, q), the Fourier-Bessel transform 
oix{s,b): 



X{s,q)= bdbjQ{qb)x{s,b). 



(5) 



In general the inputs for x(s, q) are based on analogies 
with other areas (optics, geometry,...) and/or related with 
some microscopic concepts (elementary processes), look- 
ing always for the expected connection with quarks and 
gluons, in a well established QCD bases. 

The crucial test for any input comes from the exper- 
imental data on the physical quantities that characterize 
the elastic scattering. This demands going from Eq. (5) to 
(1) and then to the differential cross section. 



^^n\Fis,q')\', 
the total cross section (Optical Theorem), 



atot = 47rlm F{s,q'^ = 0), 



the p parameter 



Re F{s,q^ = 0) 
Im F(s,g2 =0)' 



(6) 



(7) 



(8) 



experimental data and that is the point we are interested 
in. 



2.2 Eikonal zeros in the momentum-transfer space 

One aspect that exemplifies the importance of this "in- 
verse problem" is the possibility to extract model-indepen- 
dent information on the existence of eikonal zeros in the q^ 
- space. The point is that the diffraction minimum in the 
differential cross section (the dip) is generally interpreted 
as being associated with a zero in the imaginary part of 
the scattering amplitude (g^ - space) and, therefore, from 
Eqs. (1) and (5-6), it seems natural to ask if this zero in 
the amplitude is connected with, or is a consequence of, a 
zero in the imaginary part of the eikonal in the q^ - space. 
Moreover, if that is the case, what is the phenomenological 
interpretation of the zero? Since that is one of the main 
points we are interested to discuss in this work, and also 
for future reference, let us shortly review, in chronological 
order, some previous results on eikonal zeros, as well as 
some phenomenological implications (see also [T^ ). 

By the end of the sixties, the "coherent droplet model" 
introduced the idea that the eikonal in the g^-space could 
be expressed by the product of the hadronic form factors, 
which, in turn, could be assumed similar to the electro- 
magnetic form factors [Hj. Certainly, the similarity was 
thought in general geometrical terms, such as extents and 
smoothness. In this context, the imaginary part of the fac- 
torized eikonal for pp scattering reads 



and other quantities 2,4 . 

Comparisons with the corresponding data allow to check 
the ideas and, in a feedback process, to reconsider con- 
cepts and suitable phenomenological inputs, leading to 
new tests. However, to reach global and efficient descrip- 
tions with an economical number of free parameters is a 
very difficult task, even in the phenomenological context. 
One reason is because the physical quantities depend on 
both real and imaginary parts of the scattering ampli- 
tude and those, in turn, depend on the energy, momentum 
and reaction. Another reason concerns the fact that the 
model free parameters are, in general, correlated in an in- 
tricate way, turning it difficult to identify the explicit role 
or function of each parameter in the description of the ex- 
perimental data. Moreover, in most cases, one is involved 
with numerical calculation, which introduces bias, loss of 
information, and does not allow standard error propaga- 
tion. 

Certainly, the above mentioned problems can be avoided 
or, at least, minimized, if we have some kind of "empir- 
ical" or "model-independent" information on the eikonal 
directly in the q^ - space. Once having well established sta- 
tistical bases, this information may provide suitable cri- 
teria for input selections at the early stages of the model 
construction or development. This "inverse problem" con- 
cerns the determination of the eikonal from fits to the 



lmx{s,q)=Cis)Gpiq)Gpiq), 



(9) 



where C(s) is an absorption coefficient and the proton 
form factor was represented by the dipole parametrization 
for the electromagnetic form factor: 



Gp{q) = 



1 



(1+92/^2)2' 



0.71 GeV-^ 



(10) 



Therefore, in this model Im x(s, q) is positive (no change 
of sign) and the zero in the scattering amplitude has been 
interpreted as an interference effect involving both protons 
and not associated with the individual matter distribution 

m 

However, in 1975, Victor Franco presented a detailed 
fit to pp differential cross section data at -y/s = 53 GeV 
and < q^ < 5.3 GeV^ , showing that Im x may 
become negative for q'^ > 6.5 GeV^ and concluding that 
the droplet model should be modified. By that time there 
was also indication of zeros in the pion form factor |21| 
and in the proton form factor [221 ■ Moreover, in 1977, by 
means of a multipole parametrization, Maehara, Yanagida 
and Yonezawa also obtained indication of a zero in the 
imaginary part of the eikonal at q^= 6.0 GeV^, from pp 
scattering at ^/s = 53 GeV '53'. 

From a phenomenological or model point of view, the 
problem to introduce an eikonal zero has been nicely re- 
solved by Bourrely, Soffer and Wu (BSW) in 1979, by 
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means of the following parametrization for the form fac- 
tors in the - space [Tllll^ : 



3 Experimental data, Strategies and 
Ensembles 



[G{q)f fBSwiq) 



1 1 


1 


'1 - q^ja^' 






1 ^ q^jd^ _ 



(11) 



In this Section we first refer to the experimental data to 
be analyzed and recall the main problems and limitations 
related with model-independent extraction of the eikonal 
from fits to the differential cross section data. Based on 
this discussion, we introduce novel procedures and strate- 
gies through which some of these problems can be resolved 
or minimized, as demonstrated in the Sections that follow. 



where mi, m-i and a are free parameters. As referred in 
|24| . the product of two simple poles represents the "nu- 
clear form factor" and the function /bsw('3') "reflects the 
approximate proportionality between the charge density 
and the hadronic matter distribution inside a proton" . 
With the implemented "impact picture" , BSW obtained 
good descriptions of the experimental data with the zero 
fixed at « 3.81 GeV^ PO] and, more recently, at a} k, 
3.45 GeV^ 

In the eighties, by means of fits to -p-p data, Sanielevici 
and Valin obtained indication of a zero in imaginary part 
of the eikonal at q^ ~ 5.0 GeV^ j25j and latter, making use 
of the Amaldi and Schubert parametrization Furget, 
Buenerd and Valin showed that the position of the zero 
decreases from q^ ^ 8.6 GeV'^ to ~ 5 GeV^ as the energy 
increases from ^ = 23.5 GeV to 62.5 GeV 

As commented in our introduction, in 1997, Carvalho 
and Menon presented statistical evidence for eikonal zeros 
from analyses of scattering in the limited energy inter- 
val 23.5 < \/s < 62.5 GeV^. The analysis did not allow to 
infer an energy dependence, but an estimated position of 
the zero at ^ 7 ± 2 GeV^ [TS|. 

Recently, Kawasaki, Maehara and Yonezawa investi- 
gated the connections between the zeros of the scattering 
amplitude and the zeros of the eikonal in the q^ - space 
|15|. In a detailed study, the authors considered several 
parametrizations for the scattering amplitude and intro- 
duced interesting correlations between the eikonal zeros 
and a zero trajectory in the Utot versus q^ plots. One of 
the conclusions of the work is the indication of an eikonal 
zero at « 7 GeV^ (in agreement with their previous 
results |14ll23p and possibly no other zeros in the region 
below g2 « 20 GeV^. 

From the above short review we conclude that there 
are indications of an eikonal zero at ~ 7 GeV^, which, 
in the phenomenological context, can be associated with 
a zero in the hadronic form factor of the proton. How- 
ever, despite the importance of this result, the analyses 
we referred to present one or another kind of limitation, 
as we shall discuss in detail in Sec. |3| As a consequence, 
the exact position of the zero is not yet clear and, most 
importantly, the possible dependence of the position of 
the zero with the energy remains an open problem. To 
answer these and other questions (Sec. 3.2), by means of 
an improved model-independent extraction of the eikonal, 
is the aim of this work. 



3.1 Experimental data 

As mentioned before, for particles and antiparticles the pp 
and pp scattering correspond to the highest energy inter- 
val with available data. In the high energy region, >10 
GeV, data on Gtoti P and differential cross section are avail- 
able at = 13.8, 19.5, 23.5, 30.7, 44.7, 52.8 and 62.5 GeV 
for pp scattering pHll^Sll^UllSn] and at \fs = 13.8, 19.4, 31, 
53, 62, 546 and 1800 GeV for pp scattering pnUH] . More- 
over, differential cross section data from pp scattering also 
exists at Vi = 27.5 GeV and 5.5 < q^ < 14.2 GeV^ ^ 
and as discussed in the next Subsection, that set will play 
a fundamental role in our analysis. 

We have selected the differential cross section data 
above the region of Coulomb- nuclear interference, namely 
q^ > 0.01 GeV^ and the data at q'^ = (optical point) is 
determined from the corresponding values of atot and p: 



da 



0^) 



167r 



(12) 



In terms of the momentum transfer, beyond the op- 
tical point, the pp data cover the extended region 0.01 
GeV2 < g2 < 9.8 GeV^ (except for the data at 27.5 GeV); 
in contrast pp data are available only at 0.02 GeV^ < 9^ < 
4.45 GeV^. 



3.2 Parametrizations of the Differential Cross Sections 

Several authors have investigated elastic hadron scattering 
by means of parametrizations for the scattering amplitude 
and fits to the differential cross section data. The extrac- 
tion of the Profile, Eikonal and Inelastic Overlap functions 
in the 6-space and, in some special cases, the Eikonal in the 
q^-space, has led to important and novel results. For our 
purposes, some representative works are listed in Refer- 
ences [HITSIITB], EOl, EH, ESEiEll, 



and will be discussed in what follows. We first re- 
call that typical extracted results concern geometrical as- 
pects (radius, central opacity), differences between charge 
distributions and hadronic matter distributions, existence 
or not of eikonal zeros in the g^-space and, more recently, 
connections with pomerons, reggeons and nonperturbative 
QCD aspects. 

The basic input in all these analyses is the parametriza- 
tion of the scattering amplitude as a sum of exponentials 
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in r |14lll5ll2;-{| are exceptions) and fits to the differential 
cross section data. This parametrization allows analytical 
expressions for the Fourier transform of the amplitude, 
providing also analytical expressions for the quantities of 
interest in the 6-space. However, three kind of problems 
put serious limitations on the information that may be 
extracted through this procedure: 

(1) Experimental data are available only over finite 
regions of the momentum transfer (which in general is 
small, < 6 GeV^) and the Fourier transform demands 
integration from = to infinity. This means that any 
fit is biased by extrapolations, a problem that has been 
well put by R. Lombard |2HI '■ "...extrapolating the mea- 
sured differential cross section can be done in an infinite 
number of manners. Some extrapolated curves may look 
unphysical, but they can not be excluded on mathematical 
grounds." 

(2) The exponential parametrization allows analytical 
determination of the quantities in the 6-space and also the 
statistical uncertainties, by means of error propagation 
from the fit parameters. However, in this case, the trans- 
lation to the g^-space, Eq. (5), can not be analytically per- 
formed and neither the error propagation (through stan- 
dard procedures). As a consequence, the unavoidable un- 
certainties from the fit extrapolations can not, in principle, 
be taken into account. 

(3) Some analyses introduce additional constraints in 
the fit parameters so as to test some theoretical ideas or 
to obtain the dependence of some parameters as function 
of the energy. Although, in general, the number of free 
parameters is smaller than in "empirical" analyses, this 
procedure certainly introduces theoretical bias and can 
not be considered model- independent. 

All the analyses based on parametrizations of the dif- 
ferential cross sections face with one or more of the above 
mentioned problems. Only to exemplify these aspects, we 
recall that Lombard and Wilkin treated the experimen- 
tal data by generating an ensemble of fits with accept- 
able and defining an averaged value for the eikonal 
within one standard deviation. Although the extrapola- 
tions had been taken into account, the limitations in the 
interval of momentum transfer with available data did not 
allow to extract information on the eikonal at large or 
even intermediated values of the momentum transfer 32 . 
In the analysis by Amaldi and Schubert, statistical and 
scale errors from the experimental data on pp scattering 
have been propagated by means of a numerical method, 
but only to the impact parameter space |l26. . That is also 
the case in the work by Fearnley, who treated both pp 
and pp scattering (37) . In the analysis by Sanielevici and 
Valin [2S] and later by Furget, Buenerd and Valin (23, the 
eikonal has been determined in the g^-space, but without 
uncertainties from the fit parameters. At last, we recall 
that the parametrization by Amaldi and Schubert (also 
used in is strongly model dependent, since the only 
parameter depending on the energy is constrained by [231 



This embodies the geometrical scaling hypothesis, which 
is well known to be violated at the Collider energy {s/s — 
546 GeV). 



3.3 Strategies and ensembles 

Based on the above discussion, we have developed a model 
independent analysis aimed to optimize some aspects of 
this kind of approach and, simultaneously, to establish 
the confidence intervals of the extracted informations on 
well based statistical grounds. The main points are the 
following: 

(1) Compilation of the widest amount of experimental 
information presently available. 

(2) Choice of a parametrization and a fit procedure as 
model independent as possible. 

(3) Development of a statistical procedure in order to 
estimate the eikonal uncertainties in the momentum trans- 
fer space. 

(4) To perform a systematic investigation of the effect 
in the extracted eikonal, related with the existence or not 
of experimental data at large values of the momentum 
transfer. 

To this end, we first make use of the empirical ev- 
idence that at -/s ~ 20 - 60 GeV and for momentum 
transfer above g^ '--^ 3 GeV'^ , the pp differential cross sec- 
tion data are almost energy independent [411142) . This is 
illustrated in Fig. 1 for pp scattering at 19.5 GeV < ^/s < 
62.5 GeV. The data show agreement with a scattering 
amplitude parametrized by a sum of two exponentials in 
g2 and a fit through the CERN-MINUIT routine has 
given 



F{q) = 0.032e 



-1.06?'' 



0.0012e 



(13) 



a(s) 



o-tot(s) 



atotiV^ = 23.5 GeV) 



[l-zp(s)]. 



with xVDOF ^ 294/162 - 1.8. The result of the fit and 
the contributions from the above two components are dis- 
played in Fig. 1. This suggests that the data at ^/s = 27.5 
GeV, covering the region 5.5 < g^ < 14.2 GeV may be 
included in the analyses of individual sets at other ener- 
gies and our first point is to investigate how far can we 
go by including these data in different sets for pp and pp 
scattering. 

Our strategy is to consider two different ensembles of 
data, initially characterized and denoted as follows. 

Ensemble A : Original sets of data at each energy: 7 sets 
for pp scattering and 7 sets for pp scattering (Sec. 3.1). 

Ensemble B : Sets of ensemble A including in each one 
the data a.t ^/s = 27.5 GeV. 

Once selected a parametrization for the scattering am- 
plitude (to be discussed in what follows) the validity or 
not of ensemble B (that is, the compatibility or not of the 
data at 27.5 GeV with the original set) shall be checked 
by means of fits through the CERN-MINUIT routine gS] 
and standard statistical interpretation of the fit results 
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3.0 6.0 9.0 12.0 15.0 



q' (GeV') 

Fig. 1. Differential cross section data for pp scattering above 
= 3.0 GeV and fit through Eq. (13), with the contributions 
from the two exponential terms. 

We note that the addition of the data at y/s = 27.5 
GeV to pp data at 52.8 GeV has been previously explored 
by Sanielevici and Valin [221 ■ The novel aspect here is to 
perform a systematic investigation on the validity of this 
assumption in other energies and in pp scattering. 



4 Parametrization, fitting of the data and 
results 

4.1 Parametrization and fit procedure 

We consider the standard parametrization of the scatter- 
ing amplitude as a sum of exponentials in : 



Fis,q) -> ^a.e-'3•^^ 

The specific structure of the parametrization was deter- 
mined by means of the following model-independent pro- 
cedure. 

We have first started with the pp scattering at y/s — 
52.8 GeV, since the original set (ensemble A) corresponds 
to the highest interval in the momentum transfer with 
available data, namely 0.01 < < 9.8 GeV^ and the 
optical point, Eq. (12). The diffraction peak region, 0.01 
< q^ < 0.5 GeV^, is characterized by a change of the 
slope around q^ ^ 0.13 GeV^ 02] and the dominance 
of the imaginary part of the scattering amplitude, since 
p — 0.078 For these reasons we have represented 

the imaginary part of the amplitude by a sum of two ex- 
ponentials. From the visual examination of the data, we 
inferred initial values for the free parameters and then, 
they have been statistically determined by means of the 
CERN-MINUIT routine. Similar procedure has been used 
in the region of large momentum transfer (3.0 < q^ < 9.8 



GeV^), leading to the determination of two additional and 
independent exponential terms. In order to generate the 
diffraction minimum we added to the previous fixed re- 
sults an exponential term with negative sign. Looking for 
the simplest approach, with an economical number of free 
parameters, we have tested the possibility that two expo- 
nential terms could represent the real part of the ampli- 
tude. The only constraint corresponds to the definition of 
the p parameter, Eq. (8). With the "experimental" p value 
as input and using as initial values for the free parameters 
those determined previously, the fit led to a statistically 
consistent final result, x^/DOF = 323/196 = 1.65. 

From this result at ^/s = 52.8 GeV we fitted the pp 
data in a sequence of nearest energies, using as initial val- 
ues for the free parameters those previously obtained in 
each case. The same procedure has been applied to pp 
data beginning at y/s — 53 GeV and then going through 
the sequence of nearest energies. The values of the fit pa- 
rameters with ensemble A have also been used as initial 
values in the fits with ensemble B. 

With this procedure we have obtained a good repro- 
duction of the experimental data with the following par- 
ametrizations for the real and imaginary parts of the scat- 
tering amplitude: 

2 

Re F{s,q) ^ p^aje~^^'^\ (14) 



Imi^(s,g) = ^a,e-''^'', (15) 
i=i 

where 

M=^^Ea„ (16) 

cxj,l3j (j=l,2,...,n), are real free parameters and p{s) is 
the value extracted from experiments at each energy. The 
number of exponentials terms depends on the reaction and 
data set as shown in what follows. 

It should be noted that, with this procedure, the fit 
parameters are completely free, since, here, we are not 
interested in obtaining dependences on the energy and/or 
reaction. Our aim is the best statistical result, without 
any constraint in the fit parameters. 



4.2 Fitting results 

The central values of the free parameters in all the fits 
presenting statistical consistence are displayed in Tables 
1-3, together with other statistical information. The er- 
rors in the parameters (not displayed) correspond to an 
increase of the by one unity and involve variances and 
covariances which will be discussed in the next Subsection. 
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The results from pp scattering with ensembles A and 
B appear in Tables 1 and 2. We have found that the data 
at a/s = 13.8 GeV are not compatible with the data at 
27.5 GeV (ensemble B). In the case oi pp scattering none 
of the data sets are compatible with the pp data at 27.5 
GeV (ensemble B) and the results with ensemble A are 
displayed in Table 3. Therefore, in what follows, ensemble 
B (data at ^/s = 27.5 GeV added) corresponds only to 
pp scattering at 6 energies: 19.5, 23.5, 30.7, 44.7, 52.8 and 
62.5 GeV. The fit results for pp scattering together with 
the experimental data in ensembles A and B are shown in 
Figs. 2 and 3 and those for pp scattering with ensemble A 
in Fig. 4. 

Typical contributions to the differential cross section 
from the real and imaginary parts of the scattering am- 
plitude are illustrated in Figs. 5 and 6 for pp scattering 
at 52.8 GeV and 30.7 GeV, respectively. These "contri- 
butions" , generated by the fit procedure without model 
dependence (Sec. 4.1), allow to infer some interesting "em- 
pirical" results. We see that both the real and imaginary 
parts of the amplitude present only one zero (change of 
sign), the former at small values of the momentum trans- 
fer and the later at the dip position. Also, with the excep- 
tion of the dip region, which is filled by the real part, the 
imaginary part dominates at all values of the momentum 
transfer. It is also important to note that, although not 
imposed in the fit procedure, the position of the zero of 
the real part (ensemble B) is in agreement with a theorem 
recently demonstrated by A. Martin, which states that the 
real part (even amplitude) changes sign at above « 0.1 
GeV2 



4.3 Uncertainties and error propagation 

As mentioned before, we are interested in the determina- 
tion of the confidence interval of the extracted informa- 
tion, meanly related with the existence or not of experi- 
mental data at intermediate and large values of the mo- 
mentum transfer. For this reason we made use of the same 
parametrization and fit procedure for all the sets in ensem- 
ble A, even in the cases where the available data concern 
only the diffraction peak, as shown in Fig. 4. Certainly, 
this lack of information will be mirrored in the associated 
confidence interval and that is the point we are interested 
in. 

These considerations may be quantitatively exempli- 
fied as follows. In each fit, the error matrix provides the 
variances and covariances associated with each parameter 
(the numerical values do not appear in Tables 1-3 due 
to lack of space, but are available from the authors). By 
means of standard error propagation 44 the uncertainties 
in the free parameters, Aaj, A[3j, {j = 1, 2, ...) have been 
propagated to the scattering amplitude, Eqs. (14 - 16) and 
then to the differential cross section, Eq. (6), providing 
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Fig. 2. Fit results for pp scattering with ensembles A and B. 



By adding and subtracting the corresponding uncertain- 
ties we may estimate the confidence region associated with 
all the extrapolations, which cannot be excluded on statis- 
tical grounds [3S|. A typical result with ensembles A and 
B is illustrated in Fig. 7, for pp scattering at \fs= 23.5 
GeV. We see that, as expected, the effect of adding the 
experimental data at \fs— 27.5 GeV (when statistically 
justified) is to reduce drastically the uncertainty region. 
That result will be fundamental in the extraction of the 
empirical information on the eikonal, as shown in what 
follows. 



Our fit results provide useful information for stud- 
ies in the impact parameter space, namely r{s,b) and 
Ginei{s,b). However, since we are interested here only in 
the eikonal in the momentum transfer space, we postpone 
the impact parameter analysis for a future work. 



Ensemble A 



19.5 GeV 

I L-__ 




52.8 GeV 



62.5 GeV 
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Table 1. Fitting results at each energy (in GeV) for pp scattering with Ensemble A: values of the free parameters in GeV~^, 
experimental p value, maximum value of the momentum transfer in GeV^, number of experimental points (N) and per degree 
of freedom. 





13.8 


19.5 


23.5 


30.7 


44.7 


52.8 


62.5 


ai 


-5.698 X 10"'' 


-1.465 X 10"^ 


-0.2390 


-4.115 X 10"^ 


-1.090 X 10~^ 


-2.123 X 10"^ 


-4.281 X 10"^ 


Q!2 


8.775 X 10"* 


0.3155 


3.267 


3.569 


0.6257 


1.147 


2.349 


Q!3 


5.983 


4.180 


0.2299 




3.672 


3.662 


0.1802 


Q!4 


-3.670 


-3.106 






-3.062 


-3.070 




as 


5.514 


6.588 


4.657 


4.641 


7.432 


7.020 


6.326 


/3i 


173.2 


0.7601 


1.141 


0.9121 


0.6919 


0.7983 


0.9444 


/?2 


3.711 


6.085 


8.565 


8.303 


31.77 


17.39 


11.13 


/33 


2.796 


2.341 


1.285 




2.172 


2.277 


2.832 


/94 


2.425 


2.165 






2.050 


2.165 




135 


6.428 


6.083 


4.279 


4.253 


6.094 


5.739 


5.140 


p 


-0.074 


0.019 


0.02 


0.042 


0.062 


0.078 


0.095 


2 

'7m ax 


2.82 


8.15 


5.75 


5.75 


7.25 


9.75 


6.25 


N 


100 


123 


134 


173 


208 


206 


125 


xVdof 


2.03 


2.96 


1.14 


1.00 


2.13 


1.65 


1.16 



Table 2. Fitting results for pp scattering with Ensemble B: same legend as in Tabled 





19.5 


23.5 


30.7 


44.7 


52.8 


62.5 


Ql 


-9.012 X 10"^ 


-0.1935 


-4.044 X 10"^ 


-1.014 X 10"^ 


-2.683 X 10"^ 


-4.238 X 10"^ 


a2 


0.2317 


3.619 


3.474 


0.6030 


1.198 


2.023 


as 


4.192 


0.1538 




3.699 


3.652 


0.3443 


a4 


-3.101 






-3.039 


-3.078 


9.109 X 10"^ 


as 


6.650 


4.326 


4.742 


7.404 


7.002 


6.405 


ae 


-1.572 X 10~^ 


-9.176 X 10"'* 


-1.968 X 10"^ 


-1.256 X 10"* 


-1.673 X 10"^ 


-3.263 X 10"^ 


Pi 


0.6177 


0.9112 


0.9771 


0.6685 


0.9369 


1.062 




6.055 


8.178 


8.424 


32.62 


16.94 


11.96 




2.351 


0.9123 




2.227 


2.272 


3.352 


P4 


2.167 






2.095 


2.169 


3.335 


Ps 


6.094 


4.148 


4.275 


6.135 


5.704 


5.290 


Pe 


6.323 X 10"^ 


0.3558 


0.4210 


0.2113 


0.4061 


0.4671 


P 


0.019 


0.02 


0.042 


0.062 


0.078 


0.095 


2 

'7m ax 


14.2 


14.2 


14.2 


14.2 


14.2 


14.2 


N 


153 


164 


203 


238 


236 


155 


xVdof 


2.80 


1.20 


1.28 


2.13 


2.07 


1.51 



5 Eikonal in the momentum transfer space 

5.1 Semi-analytical and numerical methods 

The first steps in going to x{s^ <?) arc the determination of 
r{s,b) and then x{s,b). Denoting for short the real and 
imaginary parts by the subscripts R and /, respectively, 
and inverting Eq. (2) we have 

X.(.,6) = tan-{^.IiM_}, (18) 



//(s,5), XRi^jb) and x/(s,5). Taking into account the 
variances and covariances in the fit parameters, error prop- 
agation gives the uncertainties for all these quantities. An 
example is shown in Fig. 8, for the case of pp scattering 
at = 52.8 GeV (ensemble A). 

As shown in Sec. 2, the imaginary part of the Eikonal 
is connected with the Inelastic Overlap Function and the 
unitarity condition, Eqs. (2 - 4). It is usually named Opac- 
ity Function and that is the quantity we are interested to 
investigate. From the fit results, together with error prop- 
agation, we have found that (see Fig. 8) 



^r-f{.s,b) + [i-rR{s,b)]^i ' [i-rn{s,b)f 



X/(g,b) = ln| , , }. (19) iiy^ 



By means of the Fourier transform, Eq. (1), the param- and therefore, the Opacity Function may be approximated 
etrization (14 - 16) provide analytical expressions for FrI-s, 6),by 
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y/s: 14 19 31 52.8 62.5 546 1800 



ai 


—1.423 


—3.627 X 10 




—3.038 X 10 




—0.2038 




a2 


7.247 


7.525 




0.9131 




2.854 




as 


2.768 


0.9795 


■ 


3.471 




0.3705 




a.4, 








—3.214 








Q!8 






8.601 


7.787 


7.377 


9.953 


16.60 


Q!6 










1.659 






Pi 


2.027 


0.6502 




0.9858 




1.519 




P2 


6.429 


6.320 




151.8 




14.27 




P3 


2.487 


2.770 




2.237 




3.235 




/34 








2.213 








P6 






5.845 


5.990 


5.457 


6.470 


8.371 


P6 










17.72 






P 


0.014 


0.029 


0.065 


0.101 


0.12 


0.135 


0.14 


2 

Qmax 


2.45 


4.45 


0.85 


3.52 


0.85 


1.53 


0.626 


N 


61 


22 


23 


52 


24 


122 


47 


xVdof 


1.00 


0.57 


1.53 


1.85 


0.72 


1.01 


0.92 




q(GeV-) 

Fig. 3. Same as Fig. 2 in the region of small momentum trans- 
fer. 




q (GeV) 

Fig. 4. Fits results for pp scattering with ensemble A. 



Xi{s,b) w In 



l-rR{s,b) 



(20) 



and the uncertainty Axi determined directly from AFji 
through error propagation. 

The last step is to obtain the Fourier transform (5), 
which, due to the structure of our parametrization, can not 
be analytically performed. At this point we first make use 
of a numerical integration through the NAG routine 07] 
and the results will be presented and discussed later. Since 
the numerical integration does not allow standard error 
propagation, we have developed the following approach, 
which we shall name semi- analytical method. 

Generically, we can expand Eq. (19) in the form 



Xiis,b)^rR{s,b)+R{s,b), 



(21) 
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Fig. 5. Proton-proton scattering at = 52.8 GeV. Con- 
tributions to the differential cross section from the real and 
imaginary parts of the scattering amplitude with ensembles A 
(up) and B (down). 



where R{s, b) corresponds to the remainder of the series. 
Performing the Fourier transform we obtain 



10" 



10" 




Xiis,q) = Fi{s,q) + R{s,q). 



(22) 



0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 
q' (GeVT 

Fig. 6. Same as Fig. 5 for proton-proton scattering at ^/s 
30.7 GeV. 



and also the errors, AR(s,b), through error propagation 
from AFrIs, b). We then generate a set of numerical points 
with the corresponding errors and making use of the CERN- 
MINUIT routine this set of points with errors, R{s, b) ± 
AR{s, b), was fitted by a sum of gaussians 



Since the amplitude Fi{s,q) and errors AFi{s,q) arc di- 
rectly given by the fits, our task concerns the evaluation 
of 



%t(«'^) = E^^-^" 



(25) 



R{s,q) 



bdbJo{qb)R{s,b), 



(23) 



with the corresponding errors, AR(s,q), and this is the 
central point of the method. First, from Eqs. (20) and 
(21), the quantity R{s,b) can be evaluated 



R{s,b) = ln[ 



l-rRis,b)' 



rR{s,b), 



(24) 



A typical result is displayed in Fig. 9. With this par- 
ametrization, R{s, q) in Eq. (23) can be analytically eval- 
uated and also the errors, AR{s,q), may be estimated 
through the propagation of the errors from Aj^ Bj, as 
given by the routine. At last, Eq. (22) leads to x/(s,q) and 
the error propagation provides Axi{s,q). This method 
was used in Ref. [27j in order to determine the eikonal 
X/(s,q). The novel aspect of our approach is its use in the 
estimation of uncertainties. 
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Fig. 7. Elastic pp scattering at ^/s — 23.5 GeV. Extrapolations 
from the fits with ensembles A and B. 



10" 



10 



10"" 



10"' 



10-' 



o 9 e e , 



r,i », 



'••lli 

! 



0.0 
10' 

10° 

10"' 

10" 



1.0 2.0 
b(fm) 



3.0 



10"' 



10 



Xr ^ I I ^ 1 



0.0 



1.0 2.0 
b(fm) 



3.0 



Fig. 8. Real and imaginary parts of the Profile and Eikonal 
functions for pp at ^/s = 52.8 GeV, with ensemble A. 




b(fm) 

Fig. 9. Parametrization for the generated remainder R{s,b) 
by means of Eq. (25), from pp scattering at 62.5 GeV. 



5.2 Zeros and the eikonal at = 

In this Subsection, we concentrate on two aspects of the 
eikonal in the momentum-transfer space: the existence of 
zeros (change of sign) and the results for xi{s, 9 = 0) from 
pp and pp scattering. Here we only present the results and 
stress some aspects, postponing discussions and physical 
interpretations to the next Section. 

In order to investigate the position of the zeros and, 
mainly, to determine the uncertainties in its values, we 
follow Ref. [23 and consider the expected behavior of xi 
at large q"^, namely xi Q^^- In FigS- 10 to 14 we plot 
the quantity q^xii^il) function of for several sets 
analyzed and obtained with both the numerical and the 
semi-analytical methods. 

The results from pp scattering with ensembles A and B 
are shown in Fig. 10 in the case of the numerical method 
and in Fig. 11 with the semi-analytical method. In the 
last case, the shaded areas correspond to the uncertainties 
obtained from error propagation. 

Figure 11 shows clearly the role and the effect of data 
at large values of the momentum transfer. In fact, within 
the uncertainties, ensemble A shows evidence for the change 
of sign only at y/s = 44.7 GeV and 52.8 GeV, which corre- 
spond to the sets with the largest q^ interval with available 
data (see Fig. 3). On the other hand, with ensemble B, we 
find statistical evidence for the change of sign at all the 
energies investigated. From these plots we can determine 
the position of the zeros and the associated errors from 
the extrems of the uncertainty region (in general not sym- 
metrical). The position of the zero can also be obtained 
from the numerical method (Fig. 10), but without uncer- 
tainties. 

In the case of pp scattering the semi-analytical method 
does not provide any evidence for change of sign as illus- 
trated in Figs 12, 13 and 14. As mentioned before, the 
pp data at y/s — 27.5 GeV are not statistically consistent 
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Fig. 10. Imaginary part of the Eikonal in the momentum 

transfer space multiplied by g*, obtained by moans of the nu- 
merical method and ensembles A and B for pp scattering. The 
numbers refer to the center-of-mass energy in GeV. 
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Fig. 11. Imaginary part of the Eikonal in the momentum 

transfer space multiplied by g*, obtained by means of the semi- 
analytical method and ensembles A and B for pp scattering. 
The energies are the same as indicated in Fig. 10. 



with the pp data and, for this reaction, the regions with 
available data are very small, q"^ < 4.5 GeV^. 

We conclude that only the pp data from ensemble B 
provide statistical evidence for change of sign. The results 
obtained with both methods are shown in Fig. 15. We ob- 
serve a systematic difference in the position of the zeros as 
obtained with the numerical and the semi-analytical meth- 
ods, the former giving values about 0.9 GeV^ below those 
obtained with the later. As the energy increases from 19.5 
GeV to 62.5 GcV, the position of the zero decreases from 
8.2 to 5.7 GcV^ with the numerical method and from 8.9 
to 6.4 GeV^ with the semi-analytical method. Although 
the decreasing is not smooth, the general trend, with both 
methods, favor that behavior, indicating that the position 
of the zero decreases as the energy increases. 

As will be discussed, another quantity of interest is the 
value of the imaginary part of the eikonal at = 0. In 
this case, all the results obtained with ensembles A and B 
and through both numerical and semi-analytical methods 
are exactly the same for the central values (the numerical 



method docs not provide the uncertainties). The results 
from pp and pp scattering are shown in Fig. 16 and will 
be discussed in Sec. 6.3. Here we only note that the gen- 
eral dependence of xii^i <?) on the energy from pp and pp 
scattering, is similar to the behavior of the total cross sec- 
tions, afgf (s) and af^^ (s) . In fact that is expected since, in 
first order, Eq. (1) reads Fj{s, q) « X7(s, q) and, from the 
optical theorem, atot{s) = 47rFr(s, q = 0). We shall return 
to this point in Sec. 6.1. 



6 Discussion 

Our main novel results are displayed in Figs. 11, 15 and 

16, and refer, respectively, to the behavior of the eikonal in 
the q^ - space (semi-analytical method) , the position of the 
zero as function of the energy and the eikonal at q^ = 0. 
In this Section we discuss the applicability of these results 
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Fig. 13. Imaginary part of the Eikonal in the momentum 
Fig. 12. Imaginary part of the Eikonal in the momentum transfer space muhiphed by g^ obtained by means of the nu- 
transfer space muhiplied by g*, obtained by means of the nu- merical and the semi-analytical methods with ensemble A for 
merical and the semi-analytical methods with ensemble A for pp scattering at 53 GeV. 
pp scattering at 19.5 GeV. 



6.1 Eikonal zeros 



in the phenomenological context, in connection with those 
obtained by other authors. 

First, we must note that to find a clear dynamical 
origin for these "empirical" properties of the eikonal is 
obviously a very difficult task. For that reason, we shall 
consider here a particular framework, which, despite its 
simplicity, is suitable for the kind of points we are inter- 
ested to raise. We shall refer to the Multiple Diffraction 
Theory by Glauber the Chou-Yang model 8 and 

the impact picture by Bourrely, Soffer and Wu 10,22,24 . 
We understand that the essential ideas to be discussed 
may be extended to more realistic or general approaches 
and, in fact, as we shall show, some limited connections 
with nonperturbative QCD may also be inferred. The dis- 
cussion will be based on the above three aspects, namely 
the eikonal zeros, the eikonal at large and the eikonal 
at = 0. 



Let us discuss our results concerning the evidence of the 
eikonal zero and its dependence on the energy in the phe- 
nomenological context. We stress that we shall consider 
simple approaches, so that some important questions may 
be raised and/or discussed without the influence of "tech- 
nical" details. 

Let us consider the Multiple Diffraction Theory by 
Glauber. For the scattering between hadrons A and B 
the eikonal is given by [6ll9ll48| 

Na Nb 

i=\ j=i 

where Ga and Gb are the hadronic form factors, Na and 
Nb the number of constituents in each hadron and /y 
the individual elementary scattering amplitudes between 
constituents (parton-parton scattering amplitudes). If we 
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Fig. 14. Imaginary part of the Eikonal in the momentum 
transfer space multiplied by q^, obtained by means of the nu- 
merical and the semi-analytical methods with ensemble A for 
pp scattering at 546 GeV. 



consider, for simplicity, that the elementary amplitudes 
are all the same, denoted by /, and that NaNb = N, for 
pp scattering we have for the imaginary parts 



Xi{s,q)^NGl fi. 



(26) 



This expression indicates that, in principle, the zero in 
the imaginary part of the eikonal may be associated either 
with the form factor or with the elementary amplitude. 
Let us discuss both possibilities. 

In the phenomenological context (multiple diffraction 
models), the interpretation of the zero as associated with 
the elementary amplitude has been discussed, for example, 
in |48| . On the other hand, more recently, elementary am- 
plitudes have been determined from nonperturbative QCD 
|17II49| . by means of the Stochastic Vacuum Model |50| . 
The fundamental input is the gluon gauge-invariant two- 
point correlation functions, which eventually determines 
the structure of the elementary amplitude. Two variants 
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Fig. 15. Position of the eikonal zeros in the momentum trans- 
fer space as function of the energy, for pp elastic scattering 
in the interval 19.5 — 62.5 GeV: semi-analytical method (with 
errors) and numerical method (without errors). The lines con- 
necting the points have been drawn only to stress the general 
trend. 
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Fig. 16. Imaginary part of the eikonal at g = for pp and pp 
scattering as function of the energy. 



may be found in the literature, which can be distinguished 
by the behavior of the correlators in the region of small 
(physical) distances. From lattice QCD, in both quenched 
approximation (absence of fermions) and full QCD (dy- 
namical fermions included), the parametrized correlator 
has a divergent term, l/x*, at the origin |51| . In con- 
trast, the parametrization introduced by Dosch, Ferreira 
and Kramer is finite at the origin . It has been shown in 
that the elementary amplitudes determined with the 
lattice parametrization are characterized by a monotonic 
decrease of the amplitude with the momentum transfer, 
through positive values (no zeros) . On the other hand, the 
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finite correlator, by Dosch, Ferreira and Kramer, leads to 
an elementary amplitude which presents a zero at k, 
0.5 GeV^ (the position of the zero depends on the value 
of the gluonic correlation length) and then goes asymp- 
totically to zero through negative values ^7]. It should 
be noted that this formalism is intended for asymptotic 
energies (s — > oo) and small momentum transfer (typi- 
cally < C(l) GeV^), which put some limitations in the 
conclusions that may be inferred, as discussed in detail 
in |17ll49j . It has been claimed that the divergent term in 
the correlator is a perturbative effect that should not be 
included in nonperturbative calculations. However, since 
there is no conclusive answer to this question in the liter- 
ature, the possibility that the elementary amplitudes have 
no zeros can not be disregarded. 

Let us, therefore, consider the possibility that the eiko- 
nal zero is associated with the form factor, that is, assum- 
ing the positivity of the elementary amplitude as in the 
case of the parametrization from lattice QCD ■ In that 
case, a striking result is the dependence of the position 
of the eikonal zero on the energy. From a "pragmatic" or 
"empirical" point of view, for pp and pp hadronic elastic 
scattering, that dependence may be associated with the 
shrinkage of the diffraction peak, an effect experimentally 
verified when the energy increases in the region 23 GeV 
< y/s < 1.8 TeV In fact, as recalled before, in first or- 
der, Fj{s, q) « Xii^t 9) and, as it is known, the diffraction 
peak is dominated by the imaginary part of the ampli- 
tude. This possibility brings novel insights, since the main 
point is the implication in hadronic form factors depend- 
ing on the energy, which is an old phenomenological con- 
jecture. Despite of limited theoretical foundation, it has 
been shown, in the past, that the introduction of energy 
dependence in form factors leads to good descriptions of 
the experimental data [53.54| . In particular, pp and pp 
elastic scattering data can be well described by means of 
both geometrical models [SS] and hybrid Regge-dual mod- 
els [56 . We understand that to explore this dependence 
in well founded theoretical grounds (with analyticity and 
crossing taken explicitly into account) may lead to novel 
results in the investigation of the soft diffractive processes 
in general. 

As commented in |15| . the eikonal zero may indicate 
the existence of two components in the hadronic elastic 
process, one dominating the region of small momentum 
transfer (long range) and the other the region of large 
momentum transfer (short range) . Roughly, these compo- 
nents could be associated with the regions before and after 
the position of the zero, respectively. If that is the case the 
characteristic of the long range component is the positiv- 
ity of the eikonal in the q^ - space, in contrast with its 
negative values in the short range region. We shall discuss 
this later case in the following Subsection. 

All the previous discussion was based on high-energy 
hadronic interactions {pp elastic scattering) and, there- 
fore, the form factor concerns the hadronic structure of the 
proton. To end this subsection, let us discuss some results 
recently obtained from experiments on elastic electron- 
proton scattering and related to the electromagnetic struc- 



ture of the proton. Although, presently, discrepancies from 
two different experimental techniques characterize the re- 
sults, future experimental and theoretical developments 
may bring new insights on possible connections between 
hadronic and electromagnetic form factors, as discussed in 
what follows. 

The experiments on elastic e—p scattering provide in- 
formation on the ratio R = ^jLpG Ep{q^) I G Mpio^) , where 
Gspiq^) and GMpiq^) are the Sachs electric and magnetic 
form factors and /ip the proton magnetic moment. The tra- 
ditional method used to extract the form factors is based 
on the Rosenbluth separation technique ^7j and fits to 
data on R, as function of the momentum transfer, have 
yielded a scaling behavior: Rk, 1 |58j . 

Recently, this ratio has been measured at the Jeffer- 
son Laboratory by means of the polarization transfer tech- 
nique {ep — !■ ep) and the results indicated a decrease of 
the ratio from 0.97 to 0.27, as the momentum transfer in- 
creases from 0.5 to 5.5 GeV^ (nonscaling behavior) |59l 
I60II61| . Extrapolation from empirical fits indicates a zero 
(change of sign) at « 7.7 GeV'^ (61^ and a recent phe- 
nomenological description of these data, in the context of 
a Regge parametrization for Generalized Parton Distribu- 
tions yielded a zero in the electric form factor aX q^ k. 
GeV^ ■ Certainly, these results suggest some possible 
connections between the inferred position of the zero in 
the electric form factor and the position the eikonal zeros 
at ~ 6 - 9 GeV^ (Figure 15), which, as discussed before, 
may be associated with the hadronic form factor. 

However, new measurements performed at Jefferson 
Lab through the Rosenbluth technique have confirmed the 
scaling behavior and the same result has been ob- 
tained in more recent and improved measurements j54| : 
global analysis of the cross section data indicates i? « 1 
in the momentum-transfer region q^: 0-6 GeV^. 

As discussed in despite the inconsistency of the 
results in the region of large momentum transfer, both 
measurements (polarization transfer and improved Rosen- 
bluth) are of comparable precision and the origin of the 
discrepancy is not clear yet (see ^64. for references on pos- 
sible theoretical corrections in development). 

Certainly this enigma must be explained before any 
attempt to conclude on the existence or not of a zero in 
the electric form factor at « 7 - 8 GeV^. However, in 
case that further experimental and/or theoretical develop- 
ments might favor the polarization-transfer form factors 
|64j . it may be important to investigate possible connec- 
tions between zeros in the electric and hadronic form fac- 
tors. 

6.2 Eikonal at large q^ 

Figure 11 (and also 10) shows another interesting aspect 
of the extracted eikonal, in the region of intermediate and 
large momentum transfer. The results indicate that, after 
the zero, the eikonal reaches a minimum and then ap- 
proaches zero through negative values. That may indicate 
a second asymptotic zero, which is obviously expected as 
q^ 00. However, this approximation to zero through 
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negative values brings some new insights in model con- 
structions. 

For example, let us return to the impact picture by 
Bourrely, SofFer and Wu. In this model, the contribution 
from the Pomeron exchange is assumed to be factorized 
in s and b, with the impact parameter dependence given 
by the parametrization (11) for the form factors in the 
- space. Since, as oo, G{q) and /bsw('?) ^ — Ij 
this parametrization qualitatively reproduces the change 
of sign in the eikonal and the limit to zero through nega- 
tive values. However, as quoted before, the position of the 
zero fixed at = 3.45 GeV^ is in disagreement with 
our extracted behavior. 

It should be noted that this model overestimates the 
differential cross section data in the region of large mo- 
mentum transfer |1UII24) , a point also raised in . If this 
effect is not only a consequence of the above factorized 
Pomeron contributions, it may be associated with the be- 
havior of the eikonal after the zero, as mentioned before, 
that is, the short range contribution. In this case, it may 
be recalled that a modified BSW factor (mBSW), given 

by 




-1 I ' ' ' ' ' ' ' 1 

5 10 15 20 

q' (GeV') 

Fig. 17. The Bourrely-SofTer-Wu term, from Eq. (11), and the 
modified form, Eq. (28) , as function of the momentum transfer 
and for = 3.45 GeV^ 



/mBSw(g) = — 47^, (27) 

1 + q^/a^ 

can also generate the above behavior of the eikonal and 
leads to good descriptions of the experimental data at 
large momentum transfer. This function has been intro- 
duced in and used in both the geometrical approach 
(with = 8.2 GeV^) I^S] and a hybrid Regge-dual model 
|56|. The differences between these two ansatz are illus- 
trated in Fig. 17 for the position of the zero at — 3.45 
GeV2. 

We understand that tests with both functions, tak- 
ing into account the possibility that the parameter de- 
pends on the energy may lead to improved descriptions of 
the experimental data. We are presently investigating this 
subject. 

6.3 Eikonal at q^ = 

Figure 16 shows the imaginary part of the eikonal at q^ — 
0, as obtained with both the numerical and semi-analytical 
methods. As in the previous Subsections, we discuss here 
some simple examples concerning the applicability of these 
results in the phenomenological context. 

Let us return to the Glauber model, expressed by Eq. 
(26). Since at q^ = the hadronic form factors are nor- 
malized by Gp = 1, the dependence on the energy from 
the imaginary part of the eikonal at = may be asso- 
ciated either with the number of participants constituents 
or with the elementary amplitude, or both. Only to treat 
a toy example, let us assume N ~ Na x Nb fixed at 3 x 3 
= 9 and consider the Optical Theorem at the elementary 
level, namely (Jeiem{s) = At: fi[s,q = 0). In that case we 
can express the elementary cross section in terms of the 
imaginary part of the eikonal: 



Creiem(s) = — X/(s,g = 0). (28) 

The values of the elementary cross sections, obtained 
through this formula, with TV = 9, from the extracted 
values of the eikonal a± q^ — are displayed in Table 
4. We see that even with this toy model and for fixed 

= 9, the elementary cross sections at the ISR energy 
region, from analysis of both pp and pp scattering, is of the 
order of 5 - 6 r7i5, a reasonable estimation of the partonic 
cross sections. Parametrizations of these cross sections as 
function of the energy have been discussed in 165! • 

In this simple example, we considered the number N 
of participants in the elementary scattering as fixed. How- 
ever, in principle, this number may also depend on the 
energy and, for example, for a fixed elementary cross sec- 
tion (typically 5-6 mb) the increase in x/(s, q = 0) may 
be associated with an increase of A^. These considerations 
are aimed only to exemplify some possible uses of the ex- 
tracted eikonal. To go on with this discussion demands, 
however, a more realistic formalism. 

7 Conclusions and final remarks 

In this work we have presented the results of analytical 
fits to pp and pp differential cross section data, in a model 
independent way. As explained, we were not interested in 
the dependence of the fit parameters with the energy, but 
only in the best statistical results, in a model- independent 
context. 

By means of both a numerical and a semi-analytical 
methods, we have determined the imaginary part of the 
eikonal in the momentum-transfer space. Based on the 
confidence region of the statistical results, we conclude 
that the eikonal presents a zero and that the position of 
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Table 4. Elementary (parton-parton) cross sections from pp 
and pp scattering, Eq. (28). The errors come from the uncer- 
tainties in the extracted values of Xii^yl = 0) through the 
semi-analytical method. 



V elem- elem 

(GeV) (mb) (mb) 



13.8 


5.50±0.05 


6.12±0.02 


19.4 


5.57±0.02 


5.96±0.08 


23.5 


5.48±0.07 




30.7 


5.71±0.03 


6.05±0.02 


44.7 


5.95±0.03 




52.8 


6.07±0.02 


6.01±0.18 


62.5 


6.14±0.07 


5.99±0.26 


546.0 




9.85±0.18 


1800.0 




14.53±0.14 



the zero, roughly, decreases from 8.5 GeV^ to 6.0 GeV^ 
as y/s increases from 20 to 60 GeV. After the zero, the 
eikonal has a minimum and then goes to zero through 
negative values. 

We have presented a critical review on several aspects 
related to analytical fits to the differential cross section 
data and also discussions on the applicability of our results 
in the phenomenological context. Although limited to very 
simple models, we understand that these aspects may be 
extended to more realistic approaches. 

As discussed, in the phenomenological context, the 
positivity of the elementary amplitudes determined from 
quenched and full QCD |i49|| suggests that the eikonal zero 
might be associated with the hadronic form factor. In that 
case, the decrease in the position of the zero as the energy 
increases (Figure 15) implies in hadronic form factors de- 
pending on the energy. We understand that investigation 
of this possibility on well founded theoretical bases may 
lead to important developments in the treatment of the 
soft diffractive processes. 

We have also called the attention to the fact that the 
position of the eikonal zeros are consistent with the re- 
sults for the electromagnetic form factor obtained from 
polarization-transfer experiments. It is expected that the 
discrepancies between the Rosembluth and polarization- 
transfer form factors may be resolved in the near future 

As we have shown, only data at large values of the 
momentum transfer can provide precise answers on the 
topical question related with the eikonal zeros and their 
dependence on the energy. In that sense, it would be very 
important if the experiments at RHIC and LHC could ex- 
tend the region of momentum transfer to be investigated; 
as stated in [151: "Such experiments will give much more 
valuable information for the diffraction interaction rather 
than to go to higher energies" . 
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